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ABSTRACT 

o 
o 

^^ I We present results of an intensive two-month campaign of approximately daily 

^ I spectrophotometric monitoring of the narrow-line Seyfert 1 galaxy Akn 564 with the 

Hubble Space Telescope (HST). The fractional variability amplitude of the continuum 
■ variations between 1365-3000 A is ~ 6%, about a factor 3 less than that found in 

typical Seyfert 1 galaxies over a similar period of time. However, large amplitude, short 
JZ^ . time-scale flaring behavior is evident, with trough-to-peak flux changes of about 18% in 

X/^ . approximately 3 days. We present evidence for wavelength-dependent continuum time 

p^ I delays, with the variations at 3000 A lagging behind those at 1365 A by about 1 day. 

f^ I These delays may be interpreted as evidence for a stratified continuum reprocessing 



o 



region, possibly an accretion-disk structure. The Lya A1216 emission-line exhibits 
■^ ' flux variations of about 1% amplitude. These variations lag those at 1365A by $3 

Q^' days, and combining this with the line width yields a virial black-hole mass limit of 

Q . 1:8 X 10^ Mq. We caution, the low amplitude LyaA1216 variations may indicate the 
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bulk of the emission region is at larger radii. This scenario affects the veracity of our 
black hole mass upper limit in an uncertain manner due to the unknown nature of 
the gas velocity field. Our mass estimate is thus unreliable, however, it is consistent 
with the independent estimate M ~ 1 x 10^ Mq of Pounds et al. (2001), based on a 
fluctuation power spectrum analysis of X-ray variability in Akn 564. The black-hole 
mass and 5100 A luminosity of Akn 564 are consistent with the hypothesis that, 
relative to Seyfert 1 galaxies, NLSls have lower black hole masses and higher accretion 
rates. Other strong emission lines, e.g., CivA1549 and HeilA1640, are constrained to 
vary with amplitudes of < 5%. This low-level of emission- line variability is different 
from most Seyfert 1 galaxies, which characteristically display variations of ~ 10% on 
similar time scales. 

Subject headings: galaxies: individual (Akn 564) — galaxies: active — galaxies: 
Seyfert — Accretion discs 



1. Introduction 

Narrow- line Seyfert 1 (NLSl) galaxies were first classified on the basis of their narrow permitted 
optical emission lines, with H/3 full-width-half-maximum (FWHM) ^2000kms~^(Osterbrock &: 
Pogge 1985). They also exhibit distinctive X-ray properties (e.g., Puchnarewicz et al. 1992; Boiler 
et al. 1996; Brandt, Mathur, & Elvis 1997; Turner 1999; Leighly 1999a, b). These include a steep 
soft excess with photon index Fgoft > 3 (photon fiux Pe oc E ), a steep hard power law with 
Thard > 2.5, and rapid, short time-scale variations (e.g., ~ 30% in 1500 sec). Their UV/optical 
properties place them at one extreme of the Boroson &: Green (1992) primary eigenvector that 
has been identified in a principal component analysis; specifically, NLSl classification correlates 
with strong optical Fell and weak [Oiii] A5007 emission. NLSls are thus potentially useful in 
identifying the underlying physics that defines the primary eigenvector. 

A number of scenarios have been posited to explain NLSl properties: 

1. NLSls may have larger broad-line region (BLR) sizes compared to Seyfert 1 (SI) galaxies. 
The steeper soft excesses imply higher ionizing fluxes than SI galaxies with comparable 
luminosities. Under certain conditions this may lead to larger BLR sizes and hence smaller 
permitted line widths, assuming they reflect virialized motions about the putative black hole 
(Guilbert, Fabian, k McCray 1983; Wandel & Boiler 1998). 

2. NLSls may be low-inclination i, nearly face-on systems (Osterbrock &: Pogge 1985). The line 
widths, attributable to orbital, virialized motions in a common plane are decreased by sini, 
and the strong soft X-ray fluxes are reconcilable with accretion disk models (e.g., Madau 
1988). However, this orientation scenario has problems explaining the low [Oiii] luminosities 
often observed (e.g., Boroson & Green 1992). 



3. NLSls may have relatively low black hole masses with higher than normal accretion rates, 
compared to SI galaxies (Pounds, Done, & Osborne 1995). The smaller line widths are 
attributable to the reduced gravitational potential in which they form, and higher accretion 
rates result in luminosities that are large for their masses. 

Balmer (H/3) BLR sizes and masses have been measured for about 40 AGNs, including five 
NLSls (Wandel, Peterson, k Malkan 1999; Kaspi et al. 2000; Peterson et al. 2000; Shemmer et 
al. 2001). These results suggest that NLSls have BLR sizes comparable to those of Sis with 
similar optical luminosities. Moreover, NLSl and SI galaxies delineate a broad mass-luminosity 
relationship with the former sources populating the low mass extremum of this relationship. This 
is consistent with the hypothesis that NLSls are undermassive black hole systems with higher 
accretion rates and/or sources viewed nearly face-on. 

To distinguish between various models requires systematic measurements of the black hole 
masses M, mass accretion rates M, broad-line region sizes -RblR) and source inclinations i of NLSls 
and other Sis. These key parameters may all be measured directly or inferred through application 
of echo mapping (reverberation) techniques (Blandford & McKee 1982; Peterson 2001) that use 
the relative responses of continuum and emission-line components to constrain strongly the nature 
of the responding emission regions. The responsivity-weighted BLR size is given approximately by 
-Rblr = TC, with r the time lag, measured from cross-correlation analysis, between the continuum 
and delayed, radiatively driven emission-line variations. Combining i?BLR with the emission-line 
velocity full- width-half- maximum VfwhM; assumed to be gravitationally determined and derived 
from the variable line profile, virial masses M follow from M = /Vpwhm^blr/C, with / a factor 
of order unity that depends on the detailed emission-line gas distribution (and which may be 
determined by measuring the velocity-dependent emission- line response). The results of Wandel 
et al. (1999), Kaspi et al. (2000), and Peterson et al. (2000) are based on these reverberation 
techniques. The magnitude and wavelength-dependence of any continuum time delays, when 
combined with the spectral energy distribution, may be used to constrain the product MM under 
certain conditions (Collier et al. 1999), hence M follows given M from the associated emission-line 
reverberation measurement. The width of the continuum delay distribution at many wavelengths 
and the velocity dependent emission-line response (e.g., Welsh &: Home 1991) may be used to 
measure the source inclination. These methods for measuring M and i have not been applied 
with much success to existing datasets, since they lack sufficient signal-to-noise and monitoring 
duration. Other methods for measuring i include accretion disk fitting to Fe Ka (e.g., Nandra et 
al. 1997) and combined UV continuum and H/3 measurements (Rokaki & Boisson 1999). 

We undertook a program of coordinated multiwavelength observations of the NLSl galaxy 
Akn 564 (z = 0.0247, Huchra, Vogeley, & Geller 1999) in the summer of 2000. Akn 564 is the 
brightest known NLSl galaxy in the 2-10 keV band, with hard X-ray luminosity Lx — 10^^'^ ergs~^ 
(Turner et al. 2001). It has a steep power-law continuum (Fhard ~ 2.5), ionized Fe K line features, 
and a steep soft excess, as described by Turner et al. (2001). A 35 day ASCA observation shows 
evidence for fiux variations of ~ 36% (Turner et al. 2001, Edelson et al. 2001). Optical monitoring 



by Giannuzzo et al. (1998) reveals H/3 variations of about 8% on time scales of about 4 years, 
rather lower than observed in other SI galaxies over similar periods (e.g., about 25% for NGC 
5548). Optical spectra reveal the presence of strong [Can] and Fell emission (van Groningen 
1993). No program of UV monitoring has been previously undertaken on this source, but §3 
presents an archival lUE observation. Our program included ASCA (Turner et al. 2001 hereafter, 
Paper I), HST (this paper), optical (Shemmer et al. 2001 hereafter. Paper III), and RXTE 
(Pounds et al. 2001) observations, and represents the most comprehensive contemporaneous 
multiwavelength study of a NLSl galaxy to date. In this paper, we present the UV continuum 
and emission-line variability results obtained with HST during the period 2000 May 9 to July 8. 
These observations were primarily intended to measure the effective size of the UV BLR of a 
NLSl, and permit measurement of a virial mass M based on multiple emission lines, and thereby 
further constrain the nature of NLSl galaxies. Akn 564 also has a rich UV absorption spectrum 
(Crenshaw et al. 1999), and discussion of our absorption- line results will be deferred to a later 
paper. In §2, we describe the observations. We discuss the properties of the mean spectra in §3. 
The continuum and emission-line variability are presented in §4 and §5, respectively. Our results 
are discussed and summarized in §6 and §7, respectively. 



2. Observations 

We observed the nucleus of Akn 564 with the Space Telescope Imaging Spectrograph (STIS) 
on the Hubble Space Telescope (HST) on 46 occasions during the period 2000 May 9 to July 8. The 
first five visits were separated by intervals of ~ 5 days. Beginning on the fifth visit (2000 May 29), 
the sampling interval was decreased to ~ 1 day for the remaining monitoring period. We obtained 
low-resolution spectra during each visit with the G140L and G230L gratings, thereby providing 
full UV coverage at a spectral resolution of approximately 1.2 A over the range 1150-1730 A and 
3.2 A over the range 1570-3150 A. We used the 52" x 0".5 slit for the low-resolution spectra to 
maximize throughput and ensure accurate absolute photometry. As Akn 564 is a point source in 
the STIS spectral images, there was no significant resolution degradation by using the wide slit. 
The exposure times were 1200 and 520 seconds for each G140L and G230L spectrum, respectively, 
except on 2000 May 29, when the respective exposure times were 1434 and 816 seconds. To 
estimate the effects of the intrinsic absorption on the UV emission-line profiles, we obtained 
high-resolution echelle spectra of the nucleus on 2000 May 29 with the E140M grating at a spectral 
resolving power of A/AA « 45, 000. The E140M spectra were obtained through the 0".2 x 0".2 
aperture during four consecutive orbits to yield a combined exposure time of 10,310 sec. The 
spectra were reduced using the IDL software developed at NASA's Goddard Space Flight Center 
for the STIS Instrument Definition Team (Lindler 1998). For the low-dispersion spectra, we used 
standard point-source processing with an extraction height (perpendicular to the dispersion) of 
0".275 to obtain flux-calibrated G140L and G230L spectra as a function of wavelength. 

Small wavelength calibration uncertainties exist between spectra of a given grating. To 



remove these effects, relative wavelength calibration of spectra for each grating is required. The 
methodology we employ is analogous to that described by Korista et al. (1995) (except we do 
not subtract a continuum fit from each spectrum). We cross-correlate the spectra near the 
emission-line peaks of LyaA1216, ClvA1549, and HeiiA1640 for the G140L grating, and near 
the peaks of Cm] A1909 and MgiiA2798 for the G230L grating. Hence, we determine the nearest 
whole pixel shifts for each spectrum relative to its appropriate mean spectrum. The majority 
of the spectra (about 90%) required either a one or zero-pixel shift, and the remaining spectra 
required a two-pixel shift. The uncertainty in our relative wavelength calibration is of order ±0.5 
pixels, i.e., about 0.6 A and 1.7 A for the G140L and G230L gratings, respectively. We made no 
attempt to intercalibrate the G140L and G230L spectra, given their difference in resolution. 



3. Mean Spectra 

The observed mean and root mean square (rms) of the 46 G140L and G230L spectra are 
presented in Figures 1 and 2, respectively. The rms spectra have been corrected for the bias due 
to measurement errors. Furthermore, these spectra have been corrected for Galactic reddening 
using E{B — V) = 0.06 mag. (Schlegel et al. 1998). We note narrow line Balmer decrement and 
He II A1640/He ii A4686 ratios suggest internal reddening is potentially large for Akn 564 (Walter 
& Fink 1993; Paper HI), i.e., E(B — V) ~ 0.2 mag. Determination of the internal reddening in 
Akn 564 will be discussed in a forthcoming paper, Crenshaw et al. 2001. The mean spectra are rich 
in both emission and absorption features. The prominent emission features have been identified 
and labelled in Figs. 1 and 2. Many of the unlabelled, poorer contrast features may be associated 
with Fell emission. Intrinsic UV absorption in LyaA1216, NvA1240, SiivA1397, and CivA1549 
is present, along with Galactic absorption in lines such as SiivA1397 and MgiiA2798 (Crenshaw 
et al. 1999). We have made no attempt to identify the complex of absorption features shortward 
of LyaA1216, since this spectral region is affected by time-dependent residuals due to imperfect 
subtraction of geocoronal Lya emission. We defer discussion of the rms spectra to §5. 

In order to quantify the prominent emission-line characteristics, we determine pure emission- 
line spectra from which appropriate mean spectral measurements are then made. This was done 
through cubic spline interpolations over the absorption features and subtraction of a power-law fit 
to the continuum. Our spectral measurements (to be discussed in this section) assume zero internal 
reddening. The continuum is defined by four nominally line-free bands chosen by visual inspection 
of the mean spectrum, Fig. 1; 1155-1180A, 1350-1380A, 1460-1500A, and 1620-1660A. These 
bands only approximate the true continuum level, since Balmer recombination continuum. Fell 
blends and other weak broad emission contaminate them to some extent. The rms spectrum, Fig. 1, 
indicates these bands are not notably biased by line variability. The best-fit continuum is described 
by Fa = A;(A/1000A)° with a = -0.88 ± 0.01 and k = (1.56 ± 0.01) x lO-^'' ergss"! cm'^ A~\ 
The continuum fit uncertainties are purely statistical and should be treated with caution, since 
systematic errors due to reddening, for example, are non-negligible. This UV continuum is redder 



than Fx oc A^^'^^, typically observed in quasars (O'Brien et al. 1988). This may be due to, for 
example, disk irradiation (assuming the UV emission arises in an accretion disk) and/or internal 
reddening. We choose to extrapolate this continuum into the G230L spectral range, since here the 
continuum level is poorly defined, as the aforementioned contaminants strengthen. We note, in 
this regime, the extrapolated continuum is about 30% below the local, pseudo-continuum defined 
largely by the Fe II blends. The validity of the above continuum fit awaits detailed consideration 
of these Fe II blends and corrections for internal reddening. 

Our emission-line spectral measurements are detailed in Table 1. For each emission line in 
column (1) we give the FWHM (FW0.5) in column (2) and flux contained therein (F0.5) in column 
(3). Similarly, columns (4) and (5) give the full width at 20% maximum intensity (FW0.2), and 
flux therein (Fo.2)- Column (6) gives the line centroid Aq^2 ) defined by Fo.2- No attempt has 
been made to deconvolve blended lines, hence column (1) refers to the primary emission line only. 
In Table 1, columns (4) through (6) are empty for Ciii]A1909 and MgiiA2798 because F0.2 is 
ill-defined. Cm] A1909 is blended with Sim] A1892 and the estimated, extrapolated continuum 
near Mg 11 A2798 is significantly below the observed level. We note that all values in Table 1 are in 
the observed frame. We chose to measure line fluxes at 20% maximum intensity to minimize the 
effects of line blending, by measuring the primary emission line flux only. 

The uncertainties in the full width and flux measurements (Table 1) are dominated by 
systematic errors attributable to uncertainties in the applied absorption (where applicable) and 
reddening corrections. We assess the uncertainty associated with our absorption corrections by 
using cubic spline and linear interpolations over the LyaA1216 absorption feature. Figure 3 
presents an illustrative example of the ambiguity in the Lya A1216 emission line profile due to the 
different interpolation schemes. We find uncertainties of ~ 15% and 8% for our line width and 
flux measurements, respectively. We do not consider the systematic uncertainty in our applied 
reddening corrections, given the aforementioned possibility of large internal reddening. For 
wavelengths ^ I8OOA, systematic errors of ~ 30% due to uncertainties in our continuum fit may 
be important, but are not considered further. 

The LyaA1216, NvA1240, and CivA1549 emission lines have FWHMs of « 2000kms-\ 
compared to typical values of ~ 5000 km s~^ for SI galaxies. Their FW0.2 values are ~ 4000 km s~^. 
For HeiiA1640, the FW0.5 and FW0.2 estimates are « 1000 and 3000 kms^^, respectively (from 
the mean G140L spectrum). 

The emission-line profiles appear symmetric about their systemic wavelengths, but this 
is difficult to quantify due to contaminating emission, such as [O iii] A1664 in the red wing of 
IIeiiA1640, and absorption features affecting many of the line profiles. The centroids, Ag*^2) of the 
observed emission-line fluxes of LyaA1216, NvA1240, Siiv-^Oiv] A 1400, CivA1549, Nm] A1750, 
and MgiiA2798 are consistent with those expected based on the systemic redshift z = 0.0247, 
given the 1.2A and 3.2A spectral resolutions of the G140L and G230L gratings. The [Oi] A1336 
and Niv]A1523 emission-line flux centroids are redshifted and blueshifted by 2.5A and 1.6A, 



respectively. The [Oi]A1336 Ag™ is possibly biased by SiiA1340 emission. The HeiiA1640 
emission-line flux centroid is blueshifted by 1.3A (~ 230kms~^). The HeiiA1640 Aq^2 measured 
from the G230L mean spectrum is most likely biased by [Oiii] A1664 contamination. Figure 4 
presents the mean Lyo;A1216, CivA1549, and MgiiA2798 emission line profiles as a function of 
velocity, with the above defined continuum fit subtracted. The emission line profile amplitudes 
have been normalized to unity. The high ionization lines of Lya A1216 and C iv A1549 and the low 
ionization line of Mg II A2798 do not appear notably blueshifted or redshifted relative to systemic 
velocities. These results are discussed in §6.1. 

We find some significant differences between emission-line flux ratios in Akn 564 relative 
to those of more typical Seyfert 1 galaxies. From Table 1, measured line ratios for Akn 564 
are (Civ/Lya) = 0.27(1.09), (Ciii]/Lya)* = 0.16(0.13), (Ciii]/Civ)* = 0.65(0.12), 
(Siiv + Oiv]/Civ) = 0.48(0.07), (Siiv + Oiv]/LyQ) = 0.13(0.08), (Mgii/Lya)* = 0.33(0.22), 
and (Heii/Civ) = 0.78(0.09). The line ratios in parenthesis are those for the SI galaxy 
NGC 5548, taken from Clavel et al. (1991). Line ratios with an asterisk denote those (for 
Akn 564) estimated using the fluxes at 50% maximum intensity, compared to 20% for the others. 
This was necessary because, as mentioned previously, F0.2 is ill-defined for C III] A 1909 and 
MgiiA2798. The uncertainties in these line ratios are dominated by systematic effects, discussed 
above, and are of order 25%. This error estimate includes an ~ 15% bias due to excluding fiux 
in the line-wings. For line ratios that include Cm] A1909 and MgiiA2798, the error estimate 
may be notably underestimated on account of the continuum fit uncertainties described above. 
Our emission-line ratios show Akn 564 exhibits weaker Civ A1549 and stronger Siiv+Oiv] A 1400 
emission compared to NGC 5548. These results are consistent with earlier results of Wilkes et al. 
(1999) and Kuraszkiewicz et al. (2000). We note our Mgii/Lya ratio of 0.33 is notably larger than 
the mean value of 0.05 derived from the small sample of NLSls by Kuraszkiewicz et al. (2000). 
This is probably largely attributable to different continuum level estimates near to Mg 11 A2798. 

Finally, we checked archival lUE observations of Akn 564 to assess whether it was in a 
comparatively low or high flux state. Figure 5 presents the lUE SWP observation of 1984 Jan 
17. For comparison purposes, our mean G140L spectrum of Fig. 1 is overlaid as the thicker line. 
Both spectra have been corrected for Galactic reddening as detailed above. The continuum and 
emission-line fluxes are in qualitative agreement. In particular, the LyQA1216 and NvA1240 line 
proflles are very similar; and suggest negligible emission-line flux differences between the 1984 and 
2000 observations. Other lUE observations from 1981 are also in qualitative agreement with our 
observations. 



4. Continuum Variability Characteristics 



4.1. Light Curves 

We select five nominal continuum bands by visual inspection of the mean and rms spectra 
presented in Figs. 1 and 2: 1350-1380A, 1460-1500A, 1620~1660A, 2070-2130A, and 2960-3040A, 
henceforth referred to by their mean observed wavelengths 1365A, 1480A, 1640A, 2100A, and 
3000A, respectively. We re-iterate that each of these continuum bands overestimates the true 
continuum level on account of various contaminants (§3). These bands are not notably biased by 
emission-line variability, and thereby provide good approximations for the continuum variability 
at these wavelengths. Light curves, describing the continuum variability during the 60-day 
monitoring period, are presented in Figure 6, with the continuum bands as labelled. We note our 
observed continuum flux measurements are not corrected for reddening. 

All the continuum regions show the same qualitative behavior. The 1365A variations may be 
described by; an increase in flux of '^ 20% during the first 15 days, followed by a decrease of ~ 16% 
over the next 5 days. These variations are followed by two flares at about JD — 2450000 ~ 1697 
and 1711. The rising trough-to-peak variations are about 18% in ~ 3 days, with similar declining 
peak-to-trough variations occuring over longer time scales of several days. Both flare events appear 
to be asymmetric, the later one more so and of longer duration. Between these two events was 
a relatively quiescent period of about 5 days. Thereafter, the variations are less pronounced, at 
about the 5% level. For the longer wavelength variations, at 2100A and 3000A, the rms variations 
are of reduced amplitude with broader asymmetric flare profiles. The peaks of the 3000A flares 
are delayed, relative to those at 1365A, by 1-2 days. The contemporaneous ASCA observations 
exhibit larger amplitude correlated variations (including similar but narrower twin flare events), 
possibly delayed by about 0.4 days (Papers I and III). We note the continuum variations may be 
slightly undersampled during the intensive monitoring period of JD — 2450000 w 1694-1734, and 
are definitely undersampled for prior times, given that ~ 10% flux amplitude changes on time 
scales of 1 day are evident (see below). 

We characterize the 1365-3000 A variations by measuring two common variability parameters: 
the ratio of the maximum to minimum flux i?max and the amplitude of the intrinsic variability 
relative to the mean flux -Fvar- The latter is corrected for the measurement errors e, 

F.,r = |v'(^i-A2), (1) 

where F, a-p, and A^ = {l/N)J2i=i^^ ^'^^ respectively the standard mean flux, rms flux, and 
mean square of the measurement errors, with A^ = 46 the number of data points in the light curve 
(Edelson, Krolik, & Pike 1990, Rodriguez-Pascual et al. 1997b). The results are presented in Table 
2 for each of the light curves, as listed in column (1). Columns (2)-(5) give, respectively, the mean 
flux F, the rms flux ap, F^^.^^ and -Rmax- Both Fvar and -Rmax are potentially biased by constant 
flux components, e.g.. Fell emission. However, these effects are likely to be small. On time scales 
of about 60 days, the fractional amplitude of the intrinsic variations is essentially constant between 
1365 A and 2100 A at about 6%, with the full range of variations i?max — 1 ~ 0.31. At 3000A the 



fractional variability amplitude is about 4%, with the full range of variations i?max — 1 ~ 0.19. 
This level of intrinsic UV variability is about a factor 3 less than that found in typical SI galaxies, 
which typically display fractional flux variations i<Var ^ 18% on time scales of 60 days (Collier, 
Peterson & Home 2001). 

The 1365-1640 A continuum variations show evidence for ~ 10% flux amplitude changes 
on time scales of about 1 day, as seen in the events beginning at about JD — 2450000 ~ 1696, 
1710, and 1707, and suggest a fraction of the UV continuum emitting region must come from a 
compact region of ^ 1 light day in size. The 2100-3000 A variations show evidence for reduced 
flux amplitude changes of ~ 5% on similar time scales. For comparison, in SI galaxies far UV 
variations of about 5% rms occur on similar time scales (e.g., Korista et al. 1995; Welsh et al. 
1998). These faster, larger amplitude variations suggest NLSls (at least Akn 564) exhibit more 
variability power on short (day) time scales. We confirm this through an autocorrelation analysis. 
Figure 7 presents the Akn 564 1365A (sohd line) and NGC 7469 (a SI galaxy) 1315A (dotted 
line) autocorrelation functions (ACFs). The relative steepness of the 1365 A ACF (compared to 
the 1315A ACF) indicates Akn 564's fluctuation power density spectrum is flatter than that for 
NGC 7469, and thereby exhibits more power on short time scales. We note Pounds et al. (2001) 
have shown that the X-ray variations in Akn 564 are faster than those of typical SI galaxies. 
The full-width-half-maximum of the ACFs are 3.27 and 4.93 days for Akn 564 and NGC 7469, 
respectively, and are indicative of characteristic UV variability time scales. By assuming the mass 
ratio of the sources is determined by the variability time scale ratio of 0.66, we estimate the mass 
of Akn 564 to be M ~ 5 x 10'^ Mq, given a mass estimate for NGC 7469 of 8 x 10*^ Mq (Wandel 
et al. 1999). This mass estimate (for Akn 564) is in good agreement with that derived from our 
emission-line reverberation results, discussed in S6.2. 



4.2. Cross-Correlation Analysis 

The continuum light curves of Fig. 6 clearly display correlated variations between 1365-3000 A. 
The peaks and troughs of the flare events at JD — 2450000 ~ 1697 and 1711 occur approximately 
simultaneously; as mentioned above, the 3000A flare peaks appear delayed by about 1-2 days 
relative to those at shorter wavelengths. Differences between the light curves are also evident, for 
example, decay time scales for the two flares appear comparatively longer at wavelengths > 2100A, 
and short (day) time scale 1365A variations appear washed out at 3000A. In order to quantify the 
nature of the correlations between continuum variations we undertook a cross-correlation analysis. 
Two algorithms were employed to compute cross-correlation functions (CCFs), the interpolated 
CCF (ICCF) of GaskeU & Sparke (1986) as implemented by White k Peterson (1994), and the 
Z-transformed discrete correlation function (ZDCF) algorithm of Alexander (1997). For our data, 
the results recovered by the two algorithms are in agreement. 

Figure 8 shows the CCFs obtained by cross-correlating each of the four longer-wavelength 
continuum light curves with the 1365 A light curve. The solid line and data points with error 



- lo- 
bars show the ICCF and ZDCF CCFs, respectively, and are in good agreement. The 1365-3000 A 
variations are highly correlated as evidenced by maximum values of the cross-correlation 
coefficients rmax ~ 0.9. The 1480-2100A CCFs all peak at about zero lag, thereby suggesting 
these continuum variations occur quasi-simultaneously. The 3000A CCF peaks away from zero at 
about 0.5 days. The 2100A and 3000A CCFs appear asymmetric, compared to the 1480A and 
1640A CCFs, and suggest their responses extend over a larger range of positive delays. Figure 9 
presents similar CCFs for four optical continuum regions; 4875-4915A, 5197-5237A, 6551-6622A, 
and 6935-7004A, henceforth referred to by the observed wavelengths 4900A, 5200A, 66OOA, 
and 6900A, respectively (for consistency with Paper III). The optical data are from Paper III. 
These CCFs have been computed as described above. The optical variations are correlated with 
those at UV wavelengths with rmax ~ 0.5. The decreasing maximum correlation coefficient with 
optical wavelength is likely due to progressively less coherent time delay smeared responses and/or 
increasing noise dilution (see Table 4 in Paper III). The optical CCFs clearly peak away from zero 
at about 2 days. 

Table 3 summarizes our cross-correlation results. Column (1) lists the light curve that has 
been cross-correlated with the 1365 A continuum light curve. Columns (2) and (3) give the 
centroid of the CCF Teen, as determined by the ICCF and ZDCF algorithms. The ICCF centroid is 
calculated over all points above 0.8 times the maximum cross-correlation coefficient rmax; whereas 
the ZDCF centroid is computed from all points near the peak with a cross-correlation amplitude 
above half that of the peak. Columns 4 and 5 report the time delay Tpeak) measured from the 
peak of the CCF. Columns 6 and 7 give rmax, for the ICCF and ZDCF, respectively. Column 8 
details the FWHM of the ICCF. The reported errors on Tpeak and Teen for the ICCF are the 1-a 
uncertainties as determined by the model-independent quasi-bootstrap and flux randomization 
method of Peterson et al. (1998). We find no evidence for lags between the 1365 A and 1640 A 
light curves, as evidenced by peak and centroid lag measurements consistent with zero. The 
2100 A and 3000A variations, however, lag those at 1365 A by ~ 0.3 and 1.0 day, respectively. 
Here and hereafter, we note the centroid lags since they are a less biased estimator of the size of 
the responding emission region. These lags are greater than zero at 94% and 99%, respectively. 
The optical variations at ~ 5000 A lag behind those at 1365 A by about 2 days at no less than 
99% confidence. The optical variations at wavelengths longer than 6000 A appear to lag those at 
1365 A by a similar amount, although statistically the suggested lags are not significantly different 
from zero on account of their very low amplitude (~ 1%) variability; which is comparable to the 
measurement errors. 

We investigate these wavelength-dependent continuum lags further by binning the 46 G140L 
and G230L spectra (covering 1150-3140A) into 40 A bins, and forming light curves based on the 
total flux in each bin. We cross-correlate (using the ICCF) each of these light curves with the 
1365 A light curve. The results are presented in Figure 10. In the top panel, the histogram plot 
shows the centroid lag for each bin as a function of the bin wavelength, thus producing a "lag 
spectrum" . The corresponding values of rmax are shown in the lower panel. Again, the error bars 
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for each bin in the upper panel were determined by the method of Peterson et al. (1998). The 
sohd hne represents the best-fit function r oc (A'*' — Aq), with Aq = 1365 A, 7 = 2.4 it 0.1, and 
reduced chi-squared xt — 0-0''' for v = 50 degrees of freedom, to the lag measurements. The small 
value of x^ << 1 suggests the lag uncertainties may be overestimated. There is a clear trend of 
increasing lag with wavelength throughout the UV: the variations at 3100 A lag those at 1365 A by 
about 0.9 days. The lag measurements for A ^ 2050 A are non-zero at no less than 90% confidence. 
For wavelengths shorter than 2050 A, the lag measurements are consistent with zero delay. The 
1150-3140 A variations are well correlated with rmax ~ 0.9. In the region where the G140L 
and G230L spectra overlap, i.e., ~ 1600-1700A, the lag measurements are in good agreement. 
The lag spectrum is reasonably smooth with no clear positive deviations near the emission-line 
wavelengths, as would be the case if the emission-line response time scales were measurably larger. 
Similarly, there are no negative dips in the maximum correlation coefficient spectrum, with the 
possible exception near Mgii A2798. These results support the negligible presence of emission-line 
variability to be discussed in §5. We note the dip in the maximum correlation coefficient at about 
1700 A is due to calibration uncertainties in the G230L spectra. Similarly, the slight depression in 
r at ^ I200A is due to calibration errors in the G140L spectra. We note the two bins with the 
largest centroid lag errors are those contaminated by LyaA1216 and MgiiA2798. 

We use the optical lag measurements (Table 3) to extend the UV lag spectrum. The 
UV/optical lag spectrum is presented in Figure 11. The solid line represents the best-fit function 
r oc (A''' — Aq), with Aq = 1365 A, 7 = 1.3 ± 0.1, and Xu = 0.09 for z^ = 54 degrees of freedom. 
The dotted line represents the best-fit function to the UV data alone (as in Fig. 10), i.e., with 
7 = 2.4 lb 0.1. By including the optical data, the lag-wavelength relationship appears to flatten. 



5. Emission-Line Variability 

The largely featureless rms spectra of Figs. 1 and 2 immediately show that any emission-line 
variability is of low amplitude. Figure 1 shows LyaA1216, NvA1240, Siiv+Oiv] A 1400, 
CivA1549, and HeiiA1640 variations are present with amplitudes of < 4, 5, 6, 5, and 4%, 
respectively. Figure 2 indicates Si iii A1892, C iii] A1909, and Mg II A2798 variations with amplitudes 
of < 6, 6, and 4%, respectively. For cases where emission-line variations are present or suggested, 
the bulk of the variation occurs in the core of the line. The variations at these wavelengths, for 
emission lines contaminated by intrinsic absorption, may be a superposition of both emission and 
absorption-line variability. A discussion of the latter is deferred to a later paper. We note the 
N V A1240 variations appears to be slightly asymmetric, with a deficit of response in its blue wing. 
The variations at about 1215A are spurious, and attributable to time-dependent residuals from 
imperfect subtraction of geocoronal Lya emission. 

Emission-line light curves were constructed by direct integration, summing the flux in 
specified, absorption-free wavelength regions after subtracting a continuum defined by a linear 
fit to the nearest continuum regions on either side of the emission-line region. Our attempts 
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at validating and possibly improving our direct integration light curves through spectral fitting 
with the IRAF task SPECFIT (Kriss 1994) failed. This was because the spectral complexity of 
the emission-line profiles permitted too much freedom in the fitting process such that artifical 
emission-line variations were introduced by corresponding changes in the adopted absorption 
profile. The blanket exclusion of any absorption regions from our variable emission-line flux 
regions, defined by their respective rms spectra (Figs. 1 and 2), severely inhibited our ability to 
extract the already low-amplitude variations. We were only able to extract significant emission-line 
variations for LyaA1216. The compartively smaller fluxes and larger measurement errors for the 
other emission-lines precluding extraction of significant variations. 

Figure 12 presents the Lya A1216 light curves for the monitoring period. The data points 
with error bars describe the emission-line variations between 1240-1243 A and 1247-1250 A ("light 
curve 1"), and the dashed line those between 1240-1250 A ("light curve 2"); the latter includes 
any possible variable contribution from the intrinsic Hi absorption feature. In both cases the 
continuum is defined by a linear fit between 1155-1180 A and 1350-1380 A. The 1365A light curve, 
§4.1, is scaled and vertically shifted to fit light curve 1, and is shown by the solid line. We note 
our observed emission-line flux measurements are not corrected for reddening. 

The Lya A1216 variations resemble those of the continuum, but with greatly reduced 
amplitude and appear delayed by about a few days. In each case, the flux increases by ~ 3% 
during the first 15 days, and similar flare events occur at similar times, i.e., JD — 2450000 ~ 1697 
and 1711. The contrast of any features in the light curves is poor on account of the low- amplitude 
variability, which is comparable to the measurement uncertainties. On time scales of about 60 
days, the intrinsic line variations are ~ 1%, with the full range of variation ~ 7%. Moreover, the 
Lya A1216 light curves are clearly well correlated and exhibit similar variability patterns. These 
results suggest any intrinsic H i absorption-line variability does not notably affect the emission-line 
variations. 

Cross-correlation of the Lya A1216 and 1365A variations quantifies their apparent similarity. 
Our results are presented in Figure 13 and Table 3. Figure 13 presents LyaA1216 CCFs for both 
line light curves detailed above. The solid line and filled data points represent the ICCF and 
ZDCF CCFs, respectively, as described earlier, for light curve 1, and the dashed line and open data 
points are for light curve 2. The Lya A1216 variations are correlated with the 1365 A variations 
with a maximum cross-correlation amplitude rmax ~ 0.5. The probability of exceeding r ~ 0.5 in a 
random sample of observations (of the same number, A^ = 45, as here) drawn from an uncorrelated 
parent population is about 0.1% (Bevington & Robinson 1992), thereby the LyaA1216 variations 
are correlated with those at 1365A at about 99.9% confidence. The formal significance of this 
correlation is most likely over-estimated, since not all correlated data points are independent. 
Detailed simulations will be required in order to quantify the magnitude of this bias. By directly 
cross-correlating light curves 1 and 2, we confirm they are correlated with essentially zero time 
delay. The peak and centroid lag measurements for light curves 1 and 2 are consistent to within 
the large uncertainties. Moreover, for light curve 1 the lag measurements are consistent with zero 
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delay. Our results only permit an upper limit for the size of the Lya A1216 emission region, -RgLR 
^3 light days. We caution, the low amplitude Lya A1216 variations may indicate the bulk of the 
emission region is at larger radii, as discussed in §6.2. 



6. Discussion 

6.1. Nature of the UV Broad Line Region 

A complete study of the nature of the UV BLR of Akn 564 requires detailed photoionization 
modelling of the spectral energy distribution and is beyond the scope of this paper. Here, we offer 
a preliminary discussion of our results in the context of previous work. 

The BLR of Akn 564 is characterised by smaller gas velocity dispersions (as determined by the 
emission-line velocity full- width-half-maximum) than found in typical SI galaxies (§3), as noted 
earlier by Rodriguez-Pascual et al. (1997a), i.e., about 2000kms~^compared to 5000 km s~^. These 
authors found evidence for similar magnitude UV line widths in about a dozen NLSls; see also 
Kuraszkiewicz et al. (2000). The full widths of the low and high-ionization lines are comparable, 
see Fig. 4, and interestingly similar to those of the Niii] A1750 and Cm] A1909 intercombination 
lines. For example, the full widths of HeiiA1640, Niii] A1750, Cm] A1909, and Mgii A2798 are 
1800, 1900, 1900, and 1700 kms~^, respectively (Table 1; as measured from the G230L spectra). 

As noted in §3, the broad emission-line profiles of Akn 564 are approximately symmetric 
about their systemic wavelengths. These results, and those above, contrast with UV/optical 
observations of quasars (e.g., Gaskell 1982; Wilkes 1984; Corbin 1991) and SI galaxies, showing 
that high-ionization lines are shifted blueward of low-ionization lines, which are approximately 
at systemic redshifts, and slightly broader. Leighly (2000) presented observations of two NLSls 
(IRAS 13224-3809 and IH 0707-495) which also showed the high ionization lines were strongly 
blueshifted and much broader than the low ionization Mg ii A2798 line, similar to results for the 
NLSl galaxies 1 Zw 1 and NGC 4051 presented by Laor et al. (1997) and Peterson et al. (2000), 
respectively. Colhnge et al. (2001) find the CivA1549 hue in NGC 4051 to be blueshifted by 
about 100 kms^^. These prior results are often interpreted as evidence for BLR radial motions. 
The low ionization lines may arise in an accretion disk, with the high ionization lines emitted in 
a wind launched from this disk (e.g., Collin-Souffrin et al. 1988). Our results suggest an absence 
of outflowing high ionization BLR gas with velocities ii;230 kms^^ in Akn 564. Moreover, the 
approximate symmetry of the emission-line profiles may indicate the predominance of optically 
thin broad line emission, since it is emitted essentially isotropically and traces the full velocity 
field (Netzer 1976; Davidson 1977; Shields, Ferland, & Peterson 1995). 

The level of intrinsic line variability is notably less than observed in typical Seyfert 1 galaxies 
over comparable time scales (§5). For example, NGC 5548 (Korista et al. 1995) showed Lya A1216, 
CivA1549, and HeiiA1640 variations of ~ 13%, ~ 9%, and ~ 14%, respectively. Similarly, 
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the same lines in NGC 7469 (Wanders et al. 1997) showed variations of ~ 10%, ~ 9%, ~ 14%. 
This difference may simply be a consequence of the reduced level of assumed ionizing continuum 
variations combined with geometrical dilution due to the finite BLR size. To investigate this 
scenario we convolved the observed continuum variations with thin shell broad line regions at radii 
of 6 and 2 light days (i.e., the predicted sizes of the Lya A1216 and Heii A1640 emission regions, see 
§6.2), respectively. In these cases, we find Lya A1216 and Heii A1640 variations of order 3 and 6%, 
respectively. These results suggest this type of scenario does not completely explain the reduced 
level of emission-line variability in Akn 564. On the other hand, the line-emitting gas in Akn 564 
may be fully ionized in hydrogen, and hence optically thin to the Lyman continuum. This would 
lead to greatly reduced emission-line variability, with detectable variations occurring only for 
sufficiently large continuum variations that induce significant changes in the ionization structure 
of the emission region. Similarly, this scenario alone is likely insufficient to account for the lack of 
line variability across the observed range of ionization, as evidenced by the calculations of Ferland 
&: Mushotzky (1982). However, we note these calculations are more appropriate for SI BLRs, and 
thereby may be misrepresentative of NLSls. Some combination of both the above effects appears 
a reasonable explanation for the reduced level of emission-line variability in Akn 564. 

The unusual UV spectral properties of Akn 564 (§3), combined with stronger ionizing fluxes 
relative to other SI galaxies, suggest NLSls BLRs have a physically distinct nature from those 
found in SI galaxies and quasars. Kuraszkiewicz et al. (2000) posit NLSl BLRs characterized by 
gas cloud densities tt-h ~ 10^^ cm^^, ionization parameters U ~ 0.01, and BLR radii -RgLR ^ ^^ 
light days. On the other hand, Rodriguez-Pascual et al. (1997a) posit njj ~ lO^'^cm"^, U ~ 0.3, 
and BLR radii -Rblr ~ ^ light days; under these conditions the broad line emitting gas becomes 
fully ionized in hydrogen, and thereby optically thin to the Lyman continuum. Our emission line 
profiles and variability results hint at a predominance of optically thin broad line emitting gas in 
Akn 564. In combination with a small estimated BLR size, §5, of ^ 3 light days, these conditions 
appear more readily reconcilable with those suggested by Rodriguez-Pascual et al. (1997a), as 
noted above. 



6.2. Evidence for a Black-Hole Accretion Disk System 

The wavelength-dependent continuum time delays of §4.2 suggest a stratified continuum 
reprocessing region extending a couple of light days from the central black hole, and provide 
possible evidence for an accretion-disk structure (Collier et al. 1998). Simple irradiated accretion 
disk models with radial temperature profiles T oc R~^''^ predict r oc A^'^. Flatter temperature 
profiles predict steeper r-A relationships, and are realizable with, e.g., irradiated flared accretion- 
disk models. These models may be able to explain the steeper UV lag spectrum presented in 
§4.2. The fine details are determined by the geometry and structure of the accretion disk, relative 
prominence of viscous heat dissipation and irradiation effects, and the emission physics. Vaughan 
et al. (1999) present evidence, based on X-ray data, for a strongly ionised disk in Akn 564. Given 
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the close proximity of the UV and X-ray emitting regions (§4 and Papers I and III), our results 
support their conclusions. In the case of NGC 7469, the only other source with a clear detection 
of wavelength-dependent continuum lags ( Wanders et al. 1997, Collier et al. 1998, Kriss et al. 
2000), the UV/optical lag spectrum was consistent with r oc A^'^, and is similar to that observed 
in §4.2. On the other hand, these continuum lags may be a result of diffuse continuum emission 
from broad-line emitting clouds (Korista & Goad 2000), or some combination of both these effects. 

In §5 we estimated an upper limit of 3 light days for the size of the LyaA1216 
broad line emitting region, -RgLR- ^^ may compare this estimate with that expected on 
the basis of results obtained for other AGNs by using the Balmer emission region size 
i^gLR^lui^i'^osity relationships of Kaspi et al. (2000) and Peterson et al. (2000). The 
luminosity ALa(5100A) = AFx(5100 A)47rL>2 « 2.4 x lO^^ergs"^ adopting a mean flux 
Fa (5100 A) = 4 X 10^"*^^ ergs^"*^ cm^^ A and a distance D = cz/Hq = 99Mpc (assuming 
Hq = 75kms^^Mpc~ ). This leads to an estimate for the size of the H/3-emitting region 
-^BLR ^12 light days from the i?gLR^l^™i^osity relationship shown in Fig. 6 of Peterson et al. 
(2000). From previous monitoring programs on SI galaxies (Netzer & Peterson 1997), we can 
estimate the size of the Lya A1216 emitting region to be R-q^ ~ O.Si^g^j^. Similarly, we expect 
that the highest ionization lines, NvA1240 and Heii A1640, arise in a region of extent ~ 0.2i2gLR- 
Based on the optical luminosity of Akn 564, we thus predict -Rblr ^ ^ light days and -Rblr ^ ^ 
light days. The factor of ~ 2 difference between our measured and predicted values for -RgLR i^ 

TT /Q 

probably not significant, since the intrinsic scatter in the i?BLR^lu™i^osity relationship is of order 
a factor 5. These results give us confidence our -RgLR estimate is in the right ballpark. 

By combining our iigLR ^PPer limit of 3 light days with the emission-line velocity full-width- 
half-maximum VpwHM ~ 2500 km s^^, derived from the rms emission-line profile (and assumed to 
be gravitationally determined), we estimate a virial mass upper limit from 

M = /^FWHM^BrR/G < 8 X 10^ Mq, (2) 

where for consistency with Wandel et al. (1999) and Kaspi et al. (2000) we use / = 2>/\/2. This 
mass estimate is consistent with the independent mass estimate of M ~ 1 x 10^ Mq of Pounds et 
al. (2001), based on a fluctuation power spectrum analysis of X-ray variability in Akn 564. And 
is consistent with our mass estimate based on the variability arguments of §4.1. Furthermore, we 
note that the position of Akn 564 in the AGN mass-luminosity plane defined by Fig. 7 of Peterson 
et al. (2000), is consistent with the best-fit regression line based on the NLSl galaxies alone, see 
Figure 14. This is consistent with the hypothesis that Akn 564 harbors a comparatively small 
black hole accreting with a higher accretion rate and/or is viewed more face-on than typical SI 
galaxies; as found earlier by Pounds et al. (2001). 

We conclude with a cautionary note concerning our derived upper limits for the size of 
Akn 564's Lya A1216 BLR and black hole mass. The reduced continuum to Lya A1216 variations 
ratio of ~ 4, relative to ~ 1 for the SI galaxies NGC 5548 and NGC 7469 (references as given in 
§6.1), suggests our Lya A1216 reverberations arise from a region only about 1/4 of the size of the 
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LyaA1216 region probed by reverberation in SI galaxies; and is therefore biased low. To illustrate 
the potential importance of this bias, assume that the Lya A1216 BLR is a thin shell of radius R. 
Equal time delay intervals define equal areas on the shell, and if our Lya A1216 response from < 
3 light days is from only 1/4 of the area of the shell, we underestimate -R by a factor of 4. This 
would lead to an R-^^ upper limit of 12 light days. The bias in our black hole mass upper limit 
depends on the unknown nature of the gas velocity field, and is not readily quantifiable without 
better data. 



7. Summary 

An intensive two- month monitoring program on the NLSl galaxy Akn 564 was undertaken 
with HST during 2000 May 9 to July 8. We summarize our results as follows. 

1. The fractional variability amplitude of the UV continuum variations between 1365-3000 A 
is about 6% on time scales of 60 days. This level of variability is about a factor of three 
less than that found in SI galaxies. We find evidence for fast, large amplitude continuum 
variations, e.g., trough-to-peak flux changes of ~ 18% in about 3 days, atypical of those 
found in SI galaxies that display similar amplitude variations on longer time scales. 

2. We present evidence for wavelength-dependent continuum time delays. The continuum 
variations at 3000 A lag behind those at 1365 A by about 1 day. By combining the UV 
data with the optical data described in Paper III, we find that the variations at 5200 A lag 
behind those at 1365 A by about 2 days. These delays may be interpreted as evidence for 
a stratified continuum reprocessing region, possibly an accretion-disk structure. However, 
the delays may be a result of diffuse continuum emission from broad-line emitting clouds, or 
some combination of both these effects. 

3. The LyaA1216 emission line exhibits intrinsic rms variations of about 1% on time scales of 
60 days. These variations lag those at 1365A by ^ 3 days, and combining this with the line 
width yields a putative black hole mass limit of ^8 x 10^ Mq. This calculation assumes the 
line width is determined by the gravity of the black hole. We caution, the low amplitude 
Lya A1216 variations may indicate the bulk of the emission region is at larger radii, and 
thereby the veracity of our black hole mass limit is questionable. The possible bias in our 
black hole mass estimate is uncertain given the unknown nature of the gas velocity field. 
Despite the unreliability of our mass estimate, it is consistent with the independent estimate 
M ~ 1 X 10^ M© of Pounds et al. (2001), based on a fluctuation power spectrum analysis of 
X-ray variability in Akn 564. 

4. The root-mean-square spectra suggest other strong emission line variability, e.g., of 
CivA1549 and HeiiA1640, occurs with rms amplitudes of < 5%. The low level of NLSl 
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emission- line variability is in contrast to that found in typical SI galaxies, which display rms 
flux variations of ~ 10% on similar time scales. 
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Fig. 1. — The observed mean (top panel) and root-mean-square (rms: bottom panel) spectra of the 
46 G140L spectra. These spectra have been corrected for Galactic reddening using E{B — V) = 0.06 
mag. The mean spectrum displays many emission and absorption-line features, and the solid 
line denotes the best-fit power-law continuum, see text for further details. The continuum bands 
used in this paper are as indicated. The rms spectrum shows evidence for LyaA1216, NvA1240, 
Siiv+Oiv] A 1400, CivA1549, and Heii A1640 variations with amplitudes of < 4, 5, 6, 5, and 4%, 
respectively. 
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Fig. 2. — The observed mean (top panel) and root-niean-square (rnis: bottom panel) spectra of the 
46 G230L spectra. These spectra have been corrected for Galactic reddening using E(B — V) = 0.06 
mag. The mean spectrum displays many emission and absorption-line features. The prominent 
emission-lines are marked, including the continuum bands utilized in §4.1. The rms spectrum 
shows evidence for SiiiiA1892, Cm] A1909, and MgiiA2798 variations with amplitudes of < 6, 6, 
and 4%, respectively. 

Fig. 3. — An illustrative example of the ambiguity in the Lya A1216 emission-line profile, and hence 
spectral line measurements (for which absorption corrections have been made) of §3, due to different 
interpolation schemes. The solid line histogram presents the Lya;A1216 emission-line profile, and 
the dashed and dotted histograms detail the linear and cubic spline interpolation corrections over 
the absorption feature. 

Fig. 4. — The mean LyaA1216, CivA1549, and MgiiA2798 emission-line profiles of Figs. 1 and 2 
as a function of velocity. The emission line profile amplitudes have been normalized to unity. The 
high ionization lines of Lya A1216 and C iv A1549 and the low ionization line of Mgii A2798 do not 
appear notably blueshifted or redshifted relative to systemic velocities. 

Fig. 5. — The 1984 Jan 17 lUE SWP observation of Akn 564. For comparison purposes, our mean 
G140L spectrum of Fig. 1 is overlaid as the thicker line. These spectra have been corrected for 
Galactic reddening using E(B — V) = 0.06 mag. The continuum and emission-line fluxes are in 
qualitative agreement. 

Fig. 6. — The observed continuum light curves for the 60 day monitoring period, with mean 
observed wavelength as labelled. The amplitude of the intrinsic rms variations is about 6%, and 
they all exhibit the same qualitative behaviour. 

Fig. 7.— The Akn 564 1365A (solid hue) and NGC 7469 ( a Seyfert 1 galaxy) 1315A (dotted hue) 
autocorrelation functions (ACFs), respectively. The relative steepness of the 1365A ACF (compared 
to the 1315A ACF) indicates Akn 564's fluctuation power density spectrum is flatter than that for 
NGC 7469, and thereby exhibits more power on short time scales. The full-width-half-maximum 
of the ACFs are 3.27 and 4.93 days for Akn 564 and NGC 7469, respectively, and are indicative of 
characteristic UV variability timescales. 

Fig. 8. — CCFs for four UV continuum regions, mean observed wavelengths as labelled. The solid 
line and data points with error bars detail the ICCF and ZDCF CCFs, respectively, and are in good 
agreement. The 1480-3000 A variations are strongly correlated with those at 1365A as evidenced 
by maximum cross-correlation coefficients of rmax ~ 0.9. The 1480-2100 A variations occur quasi- 
simultaneously, since their CCFs all peak at about zero lag. The 3000A CCF peaks at about 0.5 
days. 



23 



Fig. 9. — CCFs for four optical continuum regions, nominal observed wavelengths as labelled. All 
CCFs have been computed with respect to the variations at 1365 A. The solid line and data points 
with error bars detail the ICCF and ZDCF CCFs, respectively, and are in good agreement. The 
4900-6900A variations are correlated with those at 1365A, as evidenced by maximum correlation 
coefficients of rmax ^ 0.5. Moreover, they are delayed by about 2 days with respect to those at UV 
wavelengths. 

Fig. 10. — The UV lag spectrum binned to 40 A. Top panel: The solid histogram denotes the 
centroid lags. The error bars are those derived from a model-independent quasi-bootstrap and 
flux randomization method. The solid line represents the best-fit function r oc (A''' — Aq), with 
Ao = 1365 A and 7 = 2.4 it 0.1. There is a clear trend of increasing lag with wavelength. 
Bottom panel: The solid histogram denotes the maximum correlation coefficient. The 1150-3140 A 
variations are well correlated with rmax ~ 0.9. The dip in the maximum correlation coefficient at 
about 1700 A is due to calibration uncertainties in the G230L spectra. See text for further details. 

Fig. 11. — The UV/optical lag spectrum. The optical lag measurements, denoted by the open 
symbols, are derived from contemporaneous optical monitoring data presented in Paper III. The 
solid line represent represents the best-fit function r oc (A'^ — Aq), with Aq = 1365 A and 7 = 1.3±0.1. 
The dotted line represents the best-fit function to the UV data, denoted by the solid histogram, 
detailed in Fig. 10, i.e., 7 = 2.4ib0.1. By including the optical data, the lag-wavelength relationship 
flattens. 

Fig. 12. — Lya A1216 light curves for the 60 day monitoring period. The data points with error 
bars describe the emission-line variations between 1240-1243A and 1247-1250A and the dashed line 
those between 1240-1250A. In both cases the continuum is defined by a linear fit between 1155- 
II8OA and 1350-1380A. The 1365A light curve, §4.1 and Fig. 6, is scaled and vertically shifted 
to fit light curve 1, and is shown by the solid line. The emission-line variations denoted by the 
dashed line, include any possible variable contribution from the intrinsic Hi absorption feature. 
The intrinsic Lya A1216 flux amplitude variations are about 1%, and both light curves exhibit 
similar qualitiative behaviour which in turn mimics those of the continuum regions. We note the 
dashed line variations have been normalized to have the same mean value as those denoted by the 
data points. 

Fig. 13. — CCFs for the Lya A1216 emission-line. The solid line and filled data points represent the 
ICCF and ZDCF CCFs, respectively for light curve 1, and the dashed line and open data points are 
for light curve 2 (see Fig. 12). The LyaA1216 variations are correlated with the 1365 A variations 
with a maximum cross-correlation amplitude rmax ^ 0.5. Their CCFs peak away from zero at 
about 2 days, and the differences between them are not statistically significant. 
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Fig. 14. — The reverberation-based virial mass and optical luminosity relationship for AGNs. The 
filled circles are Seyfert galaxies from Wandel et al. (1999), and the open circles are QSOs from 
Kaspi et al. (2000). The large triangles are the AGNs from the same sources whose broad lines have 
widths less than ~ 2000 kms^^. Our upper limit mass estimate for Akn 564 is shown as a filled 
diamond. The best-fit regression line based on the narrow-line objects is denoted by the dashed 
line, with M oc lO-48±o.08^ 
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Table 1. Emission-Line Characteristics 



Line 


FWg.5 


F[5.5 


FWg.2 


^0.2 


Kl (A) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Lya 


2114 


5.30 


3553 


6.51 


1245.2 


Nv 


2809 


1.98 


4725 


2.56 


1270.5 


[Oi] 


788 


0.15 


1563 


0.20 


1338.1 


Siiv+Oiv] 


3270 


0.70 


4859 


0.84 


1433.9 


Niv] 


928 


0.35 


1487 


0.44 


1521.6 


Civ 


1934 


1.33 


3469 


1.76 


1586.3 


Hen 


1195 


0.91 


2585 


1.38 


1679.6 


He 11^ 


1831 


1.04 


3875 


1.44 


1682.7 


Nm]^ 


1882 


0.56 


3930 


0.80 


1794.0 


Ciii]'^ 


1920 


0.86 








Mgii^ 


1659 


1.73 









'^In kms ^. 

In 10^"*^^ ergs s~"^ cni"^. 

'^Spectral measurements from the G230L mean spectrum. 
All other measurements from the G140L mean spectrum. 

The spectra are corrected for a Galactic reddening of 
E{B — V) = 0.06, and assume zero internal reddening. 
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Table 2. Continuum and Emission-Line Variability Characteristics 



Light Curve F ^ cr | F^ar -Rma> 
(1) (2) (3) (4) (5) 



1365 A 


7.21 


0.44 


0.061 


1.31 


1480 A 


6.99 


0.46 


0.065 


1.33 


1640 A 


6.98 


0.46 


0.064 


1.33 


2100 A 


5.82 


0.30 


0.051 


1.28 


3000 A 


6.48 


0.25 


0.037 


1.19 


La A 1216^ 


16.95 


0.30 


0.015 


1.07 


LqA1216^ 


29.10 


0.42 


0.012 


1.06 



'^Units are 10^"*^^ erg cm^^ s^^A for 
continuum fluxes and 10^^^ erg cm^^ s^^ for 
the line fluxes. All light curves have iV = 46 
data points. 

Lya A1216 variations between 
1240-1243 A and 1247-1250 A. See §5 for 
further details. 

*^LyaA1216 variations between 
1240-1250 A. See §5 for further details. 
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Table 3. Cross-Correlation Results 



Band 


'cen 


peak 


^max 


FWHM« 




ICCF 


ZDCF 


ICCF 


ZDCF 


ICCF 


ZDCF 


ICCF 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


1480A 


-o.otgl 


0.0 


O.OiO.O 


0.0±0.4 


1.0 


1.0 


3.3 


1640A 


o.ii^J 


0.2 


O.OiO.l 


0.0±0.4 


0.9 


1.0 


3.0 


2100A 


0.3tSi 


0.3 


0.310-2 


O.OIOJ 


0.9 


0.9 


3.0 


soooA 


1 n+o-4 


0.7 


o.5tl:l 


o.ol^j 


0.8 


0.7 


3.5 


4900A 
5200A 
66OOA 


1 S+O-5 
1-^-0.4 

1 S+0-6 

2-2_8:3 


1.8 
1.7 
2.3 


9 9+0.2 

1 4+1-0 


2.7±1.3 

1 7+2.5 

9 1+0.8 


0.7 
0.6 
0.4 


0.6 
0.5 
0.5 


3.6 
3.7 
3.8 


6900A 


2.6+^-? 


2.4 


^■otli 


9 7+1.3 
^•'-1.8 


0.4 


0.3 


3.7 


La A 1216^ 


0.9+?-^ 


0.4 


0.710-6 


1.0+0-6 


0.6 


0.5 


5.0 


LaXUW 


2.8+0-0 


1.9 


2.710:0 


i.olii 


0.5 


0.4 


3.6 



'^In days. 

*'LyaA1216 variations between 1240-1243A and 1247-1250A. See §5 for further 
details. 

■^Lya A1216 variations bet^ween 1240-1250A. See §5 for further details. 
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